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Abstract 

fvi ■ We consider a model proposed by Derrida and Pomeau (1986) and recently studied 

by Chatterjee and Durrett (2009); it is defined as an approximation to S. Kauffman's 
boolean networks (1969). The model starts with the choice of a random directed graph on 
n vertices; each node has r input nodes pointing at it. A discrete time threshold contact 
process is then considered on this graph: at each instant, each site has probability q of 
choosing to receive input; if it does, and if at least one of its inputs were occupied by a 
1 at the previous instant, then it is labeled with a 1; in all other cases, it is labeled with 
a 0. r and q are kept fixed and n is taken to infinity. Improving a result of Chatterjee 
and Durrett, we show that if qr > 1, then the time of persistence of the dynamics is 
exponential in n. 
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1 Introduction 

In this paper we consider a threshold contact process on a random graph. Let r G N and 
q G [0, 1]. For n G N, n > r, let V n = {1, . . . , n}. For each x G V^, choose r distinct points 
yi(x), . . . ,y r {x) in V n — {x}; this choice is made uniformly among all fa"~_ x v possibilities 
and independently for each x £ W Let E n = {(yj(x),x) : x G 1^,1 < « < r} and call 
G n = (V n ,E n ) the graph thus obtained, a random directed graph on n vertices and in-degree 
equal to r. 

Once G n is chosen, it remains fixed and we consider a discrete time Markov chain with 
state space {0, 1}^™ and initial configuration £o G {0, l} Vn , which will be deterministic for all 
our purposes. Let {Bf : x G V n , t > 1} be a family of independent Bernoulli random variables 
with parameter q; given £ t G {0, l} Vn , we put 

f 1 if Bf +1 = 1 and YLi ZMx)) > 0; 
5t+iW - I otherwise. 

When i?f = 1, we say that x receives input at time i; therefore, a vertex is set to 1 if and 
only if it receives input at that time and at least one of its input vertices yi(x), . . . , y r (x) was 
set to 1 at the previous time. Given A C V n , £^ denotes the chain with initial configuration 
£q = 1a, the indicator function of A. We write £f instead of Q . P n will denote the 
probability measure both for the choice of G n and for the family {Bf} (they are of course 
taken independently). 

This setting was proposed by Derrida and Pomeau ([2]) as an "annealed approximation" 
to the less tractable S. Kauffman's boolean networks (@j). These are cellular automata with 
rules of evolution randomly chosen, intended to model the interactions of genes in a cell. We 
refer the reader to pQ for the detailed scientific background, including the relationship between 
the threshold contact process defined above and Kauffman's original model. 

For fixed n and any initial configuration £o G {0, l} Vn , with probability one the threshold 
contact process eventually reaches the absorbing configuration in which all vertices are in state 
0. The main object of investigation both in [2] and [1] is the distribution of this random time 
as a function of n, in particular asn^ oo. Define p = p(q, r) as the survival probability for 
a branching process in which individuals have probability q of having r children and 1 — q of 
having none. In [Tj the following is proved. 

Theorem 1.1 [Chatterjee and Durrett 2009] 

(i.) If qr > 1, then for every rj > there exist c > 0, b G (0, 1) such that, as n —?■ oo, 




inf P n i ^— L > p-rj) -)-l. 

0<t<e cnb \ n J 

(ii.) If q(r — 1) > 1, then for every r] > there exists c > such that, as n — >■ oo, 




In this paper we improve this result. 



Theorem 1.2 If qr > 1, then for every n > there exists c > suc/i that, as n —■ oo ; 

inf P n ^— !■ >p-»7 -»1. 

Q<t<e cn \ n J 

To explain why this result is to be expected and, in particular, the link with the mentioned 
branching process, we introduce the time dual of the model. Fix a realization of G n = (V n , E n ) 
and {Bf : x € V n ,t > 1}, define E n as the set of directed edges obtained by inverting the 
edges of E n and G n = (V n ,E n ). Fix T > and put B^' T = B%_ t for < t < T. Given 

A C V n , define £ ' =1a and, for < t < T, 



c A > T ( \ _ / 1 ^ f° r some z, i, we have yi(z) = x, £ t ' (z) = 1 and B*' 
»t+i l x i ~ 1 o otherwise. 



When Ct (-2) = 1 and B z t ' = 1, we say that z gives birth at time t, in which case yi(z), . . . , y r (^) 
will all be in state 1 at t + 1. We have the duality equation 

{#ns/0} = {|^ >T n a / 0} 

(we abuse notation associating £ € {0, 1} 17 ™ with {x € V^ : £(x) = 1}). Since we will only 
work with the dual process, we will drop the superscript T and assume that £j is defined for 
all positive times with the evolution rule defined above. 

Now, assume that n is very large with respect to r. If g is another integer that is much 
larger than r and much smaller than n, then with high probability the set 

{z £ V n : for some k < g and z\, . . . ,Zk £ V n , we have x — >■ z\ —►...—»■ Zfc — ^ ^ in G n } 

will simply be a directed tree of degree r rooted in x. Conditioning on this event, the evolution 
of |£f | up to time g will be exactly that of the branching process mentioned before Theorem 
11.11 In addition, it is not difficult to see that, without any conditioning, |£f | is stochastically 
dominated by such a process. This remarks clarify why the model exhibits two phases in exact 
correspondence with the branching process. If the expected offspring size qr < 1, then £f dies 
out faster than the corresponding subcritical branching process, and the primal £ t n rapidly 
reaches the zero state. On the other hand, if qr > 1, the above theorem states that the system 
survives for a time that is exponentially large in n, characterizing the supercritical regime. 

In the treatment of the dynamics, our proof is basically an exact repetition of that of [1]. 
What we do different is a more careful examination of the random graph. In order to argue 
that the confinement to a finite graph takes a long time to affect the dynamics, Chattterjee 
and Durrett prove an isoperimetric inequality that states that, if m is small in relation to n, 
then with high probability there are no subsets A C V n of size m such that the "influence set" 
{yi{x) : 1 < i < r, x € V n } has much less than rm elements. We push this argument further 
and control the influence set along several generations rather than only the first one. 

2 Proof of Theorem 1.2 

In all the results and proofs in this section, we assume that qr > 1. Also, once and for all 
we fix q < q such that qr > 1, 5 < (qr — 1) A 1 and g € N such that (qr — 1 — 5)(qr) 9 ~ l > 1 + 5. 

The following lemma is proved in [lj; see Lemma 2.2 and Equation (2.14) in that paper. 



Lemma 2.1 For every r] > there exist a > 0, b G (0, 1) suc/i t/«rf ; as n —> oo, 
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We now introduce some definitions and notation. Given m € N, let 
C = {l,...,m}x{l,...r} 4 , 0<i< 5 , 

For a = (o"o, . . . , o-j), cr' = ((Tq, <r^, . . . , cr') G T m , we say a -< a' either if i < j or if i = j 
and a is less than a' in lexicographic order. With this order, we can take an increasing 
enumeration T m = {a 1 , . . . , a ^+r+-+r g )my_ Theil) T? o = r^ ><7 mj. and) for i>\,T l m = 

( a (l+r+...+r l - 1 )m+l _ _ ff (l+r+...+r')ml 

Next, we endow T m with directed edges by setting 

a — >■ a' if and only if o" = (o"o, . . . , <7j), <t' = (o"o, . . . , ctj, 0^ +1 ) for some i. 

T m is thus the disjoint union of m rooted, directed and r-regular trees, each with g generations 
above the root. {0, l} Tm will be called the space of configurations. Given vertex a G T m and 
configuration tp G {0, l} Tm , ip(cr) G {0, 1} will denote the value of tp at a. 

Let A <zV n with \A\ = m. We can enumerate A = {x±, . . . x m } in the order of the indices 
of V n . Given a = (a , . . . , <jj) G T m with i > 0, let z CT = y ai {y<yi-^ (' ' ' (y<n GO) ■■■))■ Finally, 
define 

A u = {z a ' eT m :a' ^ a}. 

We now present an algorithm to construct a configuration tp = tp(A) G {0, 1} m from A. 

for j = 1 to m set tp(a J ) = 0; 

for j = m + 1 to (1 + r + . . . + r 9 )m 

if ty((<r»)o, (^')i, • • • , (o- j )i) = 1 for some /] or \z a3 £ A ffj 

then set ^{a 3 ) = 
else set ip{a 3 ) = 1 

In words, vertices are inspected in order; the roots are all set to and the other vertices are 
set to either if one of their ancestors has already been marked with a 1 or if their image 
under the map a *— > z a has never been seen before; otherwise they are set to 1. 

Lemma 2.2 Given A C V n with \A\ = m and <r n , . . . , a %k G T m , 

r n ( ma)](o = . . . = mMo") = i ) < ( m+rm re + j;' +r3m X k 

Proof. There is no loss of generality in assuming that a %a -< a %h when a < b. We then have 
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P n ( MA)](a l ><) = 1 | mA)}(° 11 ) = ■■■ = Vl>(A)]{a 1 ^) = 1 ) < — ±^- . 

i / n _ r 



Indeed, let Q lk denote the event that none of the ancestors of a %h in T m is marked with a 1 in 
ip(A). First note that {[^(A)]^ 1 *) = 1} C © lfe , because the algorithm fills all positions above 
a 1 with O's. Next, fix a m+ i, a m+ 2, . . . , o,i k -i € V n such that 

{z° m+1 = a m+1 , . . . ,z ffik ~ 1 = o^-i} C9"n{ [^)KO = • • • = [^)](^ fc_1 ) = 1 } 

(we start at m + 1 because z a ,... ,z a are deterministic, equal to the points of A). Then, 
conditioned on {z a = a m +i, . . . , z u k = aj fc _i}, there are at least n — r possible positions 
for z a , and [?/>(j4)](,z ct k ) = 1 precisely when z a k £ A a k , a set of size less than m + rm + 
. . . + r 9 m. I 

For i < g and B C T { m , write J(B) = {V <G T m : a ->■ a' for some aeBjc T^ 1 " 1 . Again 
for i < g, a set of sets £> = (2?o> B±, . . . , B{) with i?o C Z^ and -Bj+i C J(Bj) \/j € {1, . . . , i— 1} 
will be called a family in T m . Given ^ G {0, l} Tm , a family i3 = (So, . . . , Bi) in T m is called 
^-admissible if 

(i.) U; =0 ^C{V = 0}; 

(ii.) \B \>qm; Vj € {1, . . . ,i}, |Bj| > g ■ |J(Bj-i) n {^ = 0}|. 
,6 will be called "ip-good if it is ^-admissible and 

(m.) Vj e {0,1,..., i}, | J(Bj) n {-0 = 0}| > (gr-l-<5)(gr) J m. 

Finally, we say that ^ is robust if all -^-admissible families are V'-good. The next lemma shows 
the motivation for these definitions. 

Lemma 2.3 There exist c, C > such that, for every n and A C V n , 
^n (\if\ < 0- + S)\A\ | 4>{A) is robust) < Ce _c|A| . 

Proof. We define inductively a family B = (Bq, . . . , B g _i): 

-E>o = {°" £ ^Tra g : ^ CT gi ves birth at time 0}, 
B j+1 = {a € J(Bj) n {V'(^) = 0} : z a gives birth at time j + 1}. 

The definition of Bq implies that £^ Z> {z CT : o~ £ J(f?o)}- From the construction of ip(A) we 
see that a i-> z CT is injective on ./(So) n {V'(^) = 0}, so we have \£f\ > \J(B ) n {V'(^) = 0}|. 
Iterating this argument we get 

\£f\ > \J(Bj-i) n MA) = 0}|, l<j<g. (2.1) 

From this equation with j = g, property (Hi) in the definition of ip(A)-good families and the 
choice of g in the beginning of this section, we see that the result will follow from 

P n (B is not ^(^l)-admissible | ip(A) is robust) < Ce _c|A| Vn. (2.2) 

Define the events Go = {\B \ < qm}, Gj = {\Bj\ < q- \ J(Bj^) n{ip(A) =0}|},1 <j<g-l. 
We now have {ip(A) is robust, B is not 7/>( J 4)-admissible} C U g ,Z ({^(A) is robust} n Gj), so 

P n (B is not V(-4)-admissible | ip(A) is robust) < P„(G | ip(A) is robust) + 

3-1 

^2^n( Gj | ip(A) is robust, n^Gf ) . (2.3) 



In order to bound the terms of the sum, we will need the estimate 

x e (0, 1) ==> P( Bin(k,p) < xkp ) < exp{-'j(x)kp}, (2.4) 

where j(x) = xlogx — x + 1. This follows from Markov's inequality; see Lemma 2.3.3 in [3]. 
We start noting that 

Pn(Go | ip(A) is robust) = P( Bin(m, q) < qm ) < exp{— n f(q/q)qm}. (2-5) 

Next, on the event {ip(A) is robust, ri;~ G^}, the family £> J = (Bq, . . . ,-Bj-i) is ip(A)- 
admissible, hence tp(A)-good because ip(A) is robust, so \J(Bj-\) n {ip(A) = 0}| > (qr — 
1 — 5)(qr)^ l m > (qr — 1 — S)m. For Gj to occur, less than q\J(Bj-\) fl {^(^4) = 0}| can give 
birth at time j. Using (|2.4p . we get 



P n ( Gj | i>(A) is robust, H^Gf J < exp{--f(q/q)q(qr - 1 - <5)m}. (2.6) 

Putting d23|) and (USD together back in (|2T5|) . we get (|2T2|> . I 

Lemma 2.4 There exist e, d, D > suc/i i/iai 

m < en =4> P n (EU C K : |A| = m, tp(A) is not robust) < De~ dm Vn. 

Proof. Fix n G N and icK with |^4| = m. Let dj = |{^^4) = 1} D Z£J, 1 < i < g and 
d = J^ dj. Let i < g and B = (Bq, . . . , £>j) be a ^(.A)-admissible family. We have 

I -Bo I > Qm; 

| J(Bq) fl {ip = 0}\ > qrm — d\ > qrm — d; 

\B\\ > q rm — qd±; 

\J(B\) n {ip = 0}| > (qr) m — grdi — d2 > (qr) m — qrd; 

\J(Bi) n {V = 0}| > (gr) i+1 m - (§r)<di - (gr)*" 1 ^ - • • • - ?rdi_i - 4 
> (gr)* +1 m - (qrfd. 

Suppose that d < (1 + 5)m. Then, for j € {0, . . . , i} we have 

\J(Bj) n {ip = 0}| > (<?r) J ' +1 m - (gr^d > (qrY(qr - 1 - 5)m, 

which means that B is ?p(A)-good. We have thus shown that {ip(A) is not robust} C {d > 
(1 + 5)m}. From this we get 

P n (3A C V n : \A\ = m, rp(A) is not robust) < V] P n (i/>(A) is not robust) 

A:|A|=m 

(14-TH hr 9 )m 

^ E E E Pn([^(A)](a)=lV(7€i?). 

A:|A|=m d=r(l+<5)m] £>CT m :|£>|=d 



We now bound \{D C T m : \D\ = d}\ by 2' Tm ' and use Lemma 12.21 to bound the probability; 
the above is then less than 

n ) (1 + r + • • • + r°) m 2 (^+"^)- f (1 + r + • • • + r»)m \ ^ m 
mj \ n — r J 

< (™) m c™ (-) il+S)m ( n \ (1+S)m < (c(-) 5 ^ m 

here G is a constant that only depends on r, g and 5, and whose value has changed in the last 
inequality. Now it suffices to choose e such that Ge < 1. I 

Proof of Theorem 1.1. From here on the proof continues very similarly to pQ ; we present 
it for completeness. Fix n > 0, 



• take a, b corresponding to r\ as in Lemma 12. II 

• take c,C as in Lemma 12.31 and e,d,D as in Lemma l2.4| 

• define F = C + D, f = cAd, 1= |eW 2 ) ft J ; 

• assume that n is large enough so that n b < en < I and 5n b > 1. 

Define 

Si = \a log n] + ig, 

on = ( [n b \ + 1) A [en\ , < i < I. 

Given A C V n and a G N, let TL a (A) denote the first a elements of A (with respect to the 
order of V n ). For x £ V n , define 

Co" = n ao (|* ) , o+i = n ai+1 ( gt£ ) , o < i < i, 

where £ s ]+i denotes the set of vertices that at time Sj+i descend from Q at time Sj. Finally, 
define the events 

Ff = MCf) is robust}, Gf = {|£r|>(l + <5)ICfl}, < i < I; 
H x _ 1 = {e so >n b }, Hf = H*_ 1 nF?nGf, 0<i<I; 

Fix i > and assume that Hf occurs. Then, |£q| = oq and, by the definition of G§, we 

have \£ s °' °| > (1 + 5)|Co I > ICol because of the hypothesis 5n b > 1. So we have |(f | = cci, 
and arguing similarly we get |£J| = ctj for 1 < j < i + 1. Since £J C £g. for each j, we get 

l^faiognl+jgl - Kj I = "j for each J- As a consequence, 

flf C {|f / 0} (2.7) 

since [alogn] + ( I + l)g > /. We then have 

{x £ V n : |f / 0}| > (p - V )n) >V n (Hn{\{xeV n : |CI > «o}| > (p - */)«}) , 
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the reason being that, on the event in the second probability, at least (p — n)n vertices reach 
time so with more than ao descendents, and all that do so continue having more than a.- L 
descendents at times s,, for all i < I, and in particular are alive at time /, as seen in (|2.7p . 

We know from Lemma 12. II that, as n — > oo, 



{x e V n : \i x so \ > a } > {p -r])n) -► 1. 
Using Lemmas 12.31 and 12.41 we also have 



Pn(# C ) < E E [ P ^-l n (^ ) C ) + P n(#f-1 n (Gf) c )] 
x£V„ 8=0 
I 

< Y^ E { Fn (3A C V n : \A\ = at, tp(A) is not robust) 

x£V n 1=0 

+ P„ (\S^\ < (1 + <S)|C?I | -0(C) is robust)" 
- E E[ Ce ~ CQ ' + ^ e ~ dai ] < n ( (en - ra 6 )Fe~ /LnbJ + (I - en)F e - /LmJ 

ajGVn i=0 

which tends to as n — > oo. We thus have 



lim P n 

n— ¥oo 



{x€V n :tif ^0} >(p-T7)n) =1. 



Using duality and atractiveness, we obtain 



inf P r , 

0<i<7 



\St 



V n \ 



n 



>p-rj 



Is/ 



n 



> p-n 



|{x:£?^ 



n 



completing the proof. 



>p-r)\ -»-l, 
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